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Abstract 



We present a simple covariant constituent diquark-quark model for the nu- 
cleon. The nucleon is assumed to be composed of a scalar diquark and a 
quark which interact via a quark exchange. Starting from the Bethe-Salpeter 
equation, the instantaneous approximation leads to a diquark-quark Salpeter 
equation. In the Mandelstam formalism, the electromagnetic form factors of 
the nucleon are calculated for momentum transfers up to q 2 = —3 (GeV/c) 2 . 
A remarkable description of the experimental data is obtained. Especially, 
the model gives nearly the right values for the proton and (negative) neutron 
charge radii, and a qualitative description of the magnetic form factors. 
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I. INTRODUCTION 



The possible existence of diquarks within the nucleon has been discussed for a long 
time. Undoubtly, there is much experimental evidence for diquarks or diquark correla- 
tions/clustering. For a review on that subject see e.g. [T|-§[. One main argument for di- 
quarks is the possible explanation of the negative neutron mean square charge radius [0, 
which cannot be explained in a pure SU(6) quark model. 

On the other hand, of course, the analytical and numerical treatment of a two-body system 
is much less complicated than a three-body calculation. This is especially the case if one 
treats the problem relativistically using a Faddeev-type equation. In the last years, using 
the notion of diquarks, much progress has been done in this field. Especially, the description 
of baryons in the framework of the Nambu-Jona-Lasinio model has been successful. 
There, two quarks within the baryon are treated as diquark states analogous to mesons 
(both are 3 in colour space). The interaction with the third quark then is modelled by a 
quark exchange. Lichtenberg [§] proposed this short-range interaction long ago to split the 
degeneracy of the (70, + ) and (56, + ) ground state in a quark model of a SU(6) 21-plet 
diquark and a 6-plet quark. If one assumes a diquark clustering in the nucleon, it is very 
natural to assume such an interaction. In that case, none of the identical quarks is distin- 
guished although we have diquark correlation. In I] masses and diquark form factors are 



calculated within the NJL model. Cahill et al. []T2 1 have calculated the masses of scalar and 



axial-vector diquarks by using an approximate form of the Bethe-Salpeter equation for qq 
bound states. Most of the above references conclude a probable dominant role of the scalar 
diquark compared to the axial-vector one. The parameters used in the above literature will 
serve us as a guideline. 

Any relativistic model of the nucleon should of course reproduce correctly the electro- 
magnetic form factors of the nucleon for non-zero momentum transfers. Meyer |13|] solves 



numerically the BS equation for quark-diquark bound states with a quark exchange potential 
without, however, calculating any dynamical observables. Kroll et al. @,1^] successfully 



calculate form factors of the nucleon in a diquark model using light-cone distribution am- 
plitudes. Of course, their model is applicable only at higher momentum transfers. 

In our framework, form factors can be calculated in a straightforward way once the 
(covariant) amplitudes are known. We present a relativistic constituent quark model of the 
nucleon where the nucleon is assumed to be a bound state of a scalar diquark and a quark. 
Following the motivation given above, the interaction between the two particles is modelled 
by the exchange of a quark. Since we are only interested in the (ground state) nucleons, we do 
not incorporate a phenomeno logical confinement potential. We just want to show, that, by 
using only a (short-range) quark exchange potential, it is possible to describe the form factors 
with physically acceptable parameters. The fundamental equation is the Bethe-Salpeter 
equation for a bound state of a scalar particle and a fermion. Assuming an instantaneous 
interaction kernel we derive a Salpeter-type equation. This can be cast in a Schrodinger-type 
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equation. This eigenvalue equation is solved numerically by standard techniques. The formal 
covariance of the Salpeter equation and the Salpeter amplitudes is discussed extensively in 
fTEfl . We then use a straightforward method to calculate the electromagnetic currents in the 
Mandelstam formalism, analogous to |TO where meson form factors are computed. Since 



form factors are directly related to wave functions, the calculation of them is a crucial test 
for any model. 

This paper is structured as follows. In Sect.|| we derive the Salpeter equation for a 
bound state of a scalar diquark and a quark, including the normalization of the amplitudes. 
The assumed interaction between diquark and quark is outlined in Sect.|TB[ In Sect . [II C| the 



procedure solving the Salpeter equation is explained. The computation of the electromag- 
netic current and of the resulting form factors of the nucleon is discussed in Sect.|T|. The 
results for the form factors are presented in SectJIV] and compared to experimental data. 
Finally, in Sect.[V], a summary is given, and an outlook to an extension of the model and to 
possible applications. 



II. THE MODEL 

A. The Bethe-Salpeter equation 

A relativistic bound state of a scalar particle and a fermion with four momentum P is 
described by a Bethe-Salpeter amplitude 

Xp(xi,x 2 ) = (0\T(f>(x 1 )^(x 2 )\P) . (1) 
In momentum space, \p fulfills the following Bethe-Salpeter equation [117 1 



XP (p) = Af (pOSffo) / 0^(-iK(P,p,p')) XP (p') (2) 

with Af (pi) = if (pl—mi+ie) and S^Po) = i/ ('yp 2 —m 2 +ie) the usual Feynman propagators 
for a scalar particle and a quark, respectively, and — iK the irreducible interaction kernel. 
The masses in the above propagators are constituent masses of the order of a few hundred 
MeV. For the momenta we set 

p 1 = rj 1 P + p , p 2 = rj 2 P -p ; P = Pi + p 2 , P = V2P1 - V1P2 

with rji + r]2 = 1 . (3) 

The normalization condition is obtained as usual by considering the pole contribution of the 
two-particle (4-point-) Green's function (see e.g. 

G{x x ,x 2 , x[, 4) = ( I T0(x!Mx 2 )0ViM4) J ) (4) 
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to the bound state: 



ujp = \J P 2 + M 2 is the on-shell energy of the bound state. This yields 

J^fl dAp l (^^W^^( P ^^')]po^Xp(p) = 2u P (6) 

with 

H(P,p,p') = A F (p 1 )- 1 S[(p 2 )- l (2^S\p-p')+iK(P,p,p') . (7) 

In covariant form, Eq.|6] reads: 

d 4 p' _ . . 9 



J^y! dAp ! j^p(p) p ^ mp ^ ] ^ xM = 2M2 (8) 

Following Salpeter ]TJ| we neglect the time (i.e. energy) dependence of the interaction 
kernel by assuming K(P,p,p') = V(p,p') (instantaneous interaction) in the rest frame of 
the bound state. Then, one can easily perform the p° integration: Define in the cm. frame 
of the bound state 



HP)=( I ^XP(P°,P)) ■ (9) 



P=(M,0) 



Then, from Eq.^| one gets with standard techniques 

2loi \M — uoi — uo 2 M + oo\ + uj 2 ) J (lir)* 

with A ± (p) = ^MM. the usual Dirac projectors. H(p) = 7°(7P + m) is the Dirac Hamilto- 
nian. Define now 

W(p,p') = V(p,p') 1 ° . (11) 

Then, we can rewrite EqJTI] in a Schrodinger-type equation: 
(HV)(p) = MV(p) 



''-H 2 {p)m{p) + — f ^W(p,p')W') 



u 2 ^' ^ 2wi y (2tt) 3 

(r+v)*(p) . (12) 
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This is the fundamental equation of our model. H 2 (p) is again the Dirac-Hamiltonian of 
the quark. If \I> were a free Dirac spinor, = u(p), the kinetic energy term in Eq.[12] is 
just the sum of the kinetic energies of the two particles: T = uoi + uo 2 . In a similar way, one 
obtains the normalization condition for the instantaneous approximation from Eq.[5]: 

7^4 / d 'P Xp{p)P^^ F 1 {Pi)- 1 S^P2)- 1 )xp{p') = 2M 2 . (13) 
Performing the p° integration, one ends up with a familiar normalization of the functions 

/ ^(2a; 1 )$(p) 7 °$(p) = 2M . (14) 

It is easy to see that $(p) = &(p)j°. So apart from the factor 2ui, we obtain the usual 
normalization of Dirac spinors. In particular, the functions \& solving Eq.[T2] have non- 
negative norm for M > 0. We thus define a scalar product as follows (compare in the 
case of two fermions) : 

d 3 p 



(2vr); 



■(2a; 1 )0(p) 7 V(p) • (15) 



In order to obtain a pseudo-Hamiltonian 7i which is hermitian with respect to the above 
scalar product, the kernel W(p,p') has to fulfill the following condition: 

w(p,py = w(p',p) 

or V(p,p'^ = 7 ( V(p ,/ ,p)7° . (16) 

In general, taking into account also anti-particle states, the solutions of Eq.|T2| come in pairs. 
This can be seen as follows. Under charge conjugation, the (Bethe-)Salpeter amplitudes 
transform like (see Eq. [A5[ ) 

x b ^ f2 (p) = Sc^(-pT 

^% lh (fl = -S c W bij2 (-p) (17) 

with Sc = £ 7 2 . With b\, f 2 and bi, f 2 we denote the quantum numbers of the diquark and 
quark and their antiparticles, respectively. With ScH 2 {p)Sc = — t H2(—p) and the imposed 
condition 

W(p,p') = -Sc'Wi-p, -p')S c (18) 

Eq.[L2] results in 

(w Wi )(p) = -« w >1 
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B. The interaction kernel 



The interaction between the scalar diquark and the quark is assumed to be a quark 
exchange, see Fig. [1|. The diquark couples pointlike to the two quarks according to the 
Lagrangian 

C int = -ig s ^ c 7 5 r 2 V • (20) 
The kernel K(P,p,p') then is 

K(P, P ,p')xp{p') = -ig\ q S F q (p + p') X p{p') 

^V{pJ') = -g 2 dqq \{-%p + p')+m q ) , (21) 



with Lo q = J{p + p') 2 + m q the energy of the exchanged quark; set m q = m 2 (mass of the 
quark). In the instantaneous approximation we neglect the p° dependence of the quark 
propagator. g dqq is a dimensionless diquark-quark coupling constant. Obviously, V(p,p') 



fulfills the relation [16]. In addition, it does not mix states with different parities since 

V(p,p , )= 1 °V{-p,-p') 1 ° (22) 

and Eq.|TB| holds for the above W(p,p') = V(p, p')j°. 



C. Solving the instantaneous diquark-quark equation 

Since nucleons have positive parity and spin | we choose as basis in Dirac space the 
following states (s = ±|): 



e 1 M = ^/W + M 1 J 



e 2 s (p) = v% + M 12 ( _J_ ) , (23) 



with Myi = and W = yp 2 + M 12 . This ansatz neglects a relative angular 

momentum between the quark and the diquark. Both states are multiplied by inde- 
pendent radial functions which are chosen to be eigenstates of the harmonic oscillator: 
T^ip) — (-ll'^oE^o^ p 2/x exp(— \a 2 p 2 ). a is the oscillator parameter. Eq.[l2| is then 



diagonalized. The lowest positive eigenvalue is interpreted as the nucleon mass. It is ob- 
tained according to the Ritz variational principle by choosing the oscillator parameter a in 
the minimum of the curve M(a), see Fig|| The coupling g dqq is determined by the condition 



5 



M = 939 MeV/c 2 . We observe that for a pointlike coupling of the two quarks to the diquark, 
there is no stable solution. Only if one introduces a diquark form factor the minimum of 
the M(a) curve converges with increasing basis states, see Fig.[^. We choose a form factor 
of the following form: 

F{q) = exp(-A 2 g 2 ) (24) 



and put A = 0.18 fm, see Sect. IV, where we study the dependence of our results on A. 



III. ELECTROMAGNETIC FORM FACTORS OF THE NUCLEON 



In this section we calculate the electromagnetic form factors of the nucleon in the Man- 
delstam formalism jn],[23]. In this framework, any current matrix elements between bound 
states can be calculated in a covariant way In our model, the electromagnetic current 
coupling to the nucleon is simply the sum of the diquark- and quark current, see Fig.^|. For 
the coupling to the quark one obtains (put 771 = 772 = 5) 



P's' 



3 r rk (x = o) 



Ps 



J^yX P \P+ ^)lAAl^+P) XP{P) 



(25) 



Now introduce the vertex functions as amputated Salpeter amplitudes: 

T P (p) = Af(p 1 )- 1 S^p 2 r 1 XP (p) 
T P (p) = x P (p)A^p 1 )- 1 S^p 2 )- 1 . 



(26) 



From Eq.|2] follows in the instantaneous approximation 

d 3 p' 



r(p) :=T P (p)\ 



P=(M,0) 



(27T) 



v(p,p'Mp') 



and T(p) = — T'(jp) 
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(27) 
(28) 



Then, the rhs. of Eq.^ becomes 



d A p — 



P 



P 



P 



(27T) 



r P >( P + i)s> (- - v - q^s* (- - p)v P {p) a* (- + p) 



(29) 



Analogously, for the current coupling to the scalar diquark, one obtains 



P's' 
d A p 



j diquark j 



: (0) I Ps) = e d J J^XP>(p-~)Si(^ -P)-\P1+P' 1 ),XP{P) 



e d I —T P ip-l)SF(^-p)T P (p)A(C-+p)(P + 2p- q ) tl A((-+p-q) , (30) 



P 
~2 



P 
2" 
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where we assumed a coupling of the photon to a pointlike scalar particle. Since the vertex 
functions T{p) are only known in the nucleon rest frame (Eq.p7|), the vertex function of the 
outgoing nucleon has to be boosted according to the value of q 2 , see Eq.|M]. This covariant 
treatment of the vertex functions is very important for high q 2 . The p° integration in 



Eqs.[Z§pU| has the generic form 



dp 



f(p° 



n i = 1 , 6 (P -Pi±*e) 



(31) 



with f(p°) a regular function. The six poles arise from the denominators of the three 
propagators. The integral can be split in a residue part and a principal value integral. To 
obtain a hermitian current the latter should vanish. The <p p integration is trivial, since p and 
p' = p ± | have the same 4> p . The numerical treatment is further simplified by expanding 
the vertex functions in the basis of Eq.|2"5|: 



r.(p) = £ a](p)el + a 2 (p)el 



(32) 



where J2i sums over the radial basis states. 

We recall that the normalization condition (Eq.[TJ|) equals the zero-component of the currents 
at q 2 



0: 

P s 



jquark 



P S 



1 



P s 



Jo 



diquark 



P S 



2M 



(33) 



d? P i fr(p)At(-p)-f°T(p) , i» 2 -(-p) 7 °r( 



(27r) 3 2wi V {M-lox -Lu 2 y 



(M + lu 1 +u 2 ) 2 



The nucleon electromagnetic current can be written in the usual form 23 



(PV \j^0)\ Ps)=eu s ,{P') 



7, F^q 2 ) + 



2M 



U S (P) 



(34) 



with Fi and F 2 N being independent form factors. At q 2 = we have 



l{Ps 

Ff (0) = — 

e 



■quark 

Jo 



Ps 



Ps 



jdiquark 



P S 



2 M 



(35) 



which thus serves as a check of the numerical evaluation of the currents at zero momentum 
transfer, i.e. of the normalization. Usually, one defines the following form factors: 



G N E {q 2 ) = F^q 2 ) + ^F 2 N (q 2 ) 
GUQ 2 ) = F^q 2 ) + n N F 2 N (q 2 ) 



(36) 



The direction of the (outgoing) photon momentum is chosen parallel to the z-axis, and the 
initial nucleon is put in its rest frame. Then, G 1 ^ and G 1 ^ can be calculated by solving the 
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set of equations for j and j + = + Of course, the continuity equation has to be 
fulfilled, which serves as a check of the calculation: 



= d°j + d 3 j 3 = 



(37) 



In our framework the continuity equation follows directly from (/i = 0, 3) 
( P's' 1 3 , | Ps ) = (-l) 1 —' (P-s | j, | P'-s') 



1 



,1-s-s' 



(Ps \ 3, \P's') 



(38) 



where the transformation properties of Eqs. |29|j30| under time reversal and parity transfor- 
mation (App.[A]) have been used. Indeed, we find that the current is conserved numerically 
for any momentum transfers. 



IV. RESULTS AND DISCUSSION 



The parameter values used in our model are given in Tab.|. The parameter influencing 
the shape of the form factors mostly is the diquark parameter A (Eq.|24|). After determining 
the other parameters, A is chosen to give a best quantitative description of the proton electric 
form factor. The obtained value is in agreement with the commonly used diquark size in 
the range 0.1 — 0.3 fm []iyi5,|3"l. The quark and diquark masses are free parameters: the 
constituent quark mass is chosen to be in a range adopted by most quark models, see e.g. 
23 1 . The spectator diquark mass is a little less than the sum of two quarks Jl||],||]. These 



mass values agree with a binding energy of the nucleon proposed by the work of Ishii et al. Q 
of 50 to 150 MeV. Thus, the nucleon appears as a loosely bound state of quarks. As outlined 
in Sect. |II Cl the oscillator parameter is given by the minimum of the M(a) curve, and the 
coupling gd q q is adjusted by fitting the minimum to the nucleon mass of 939 MeV. Note that 
the so obtained oscillator parameter a = 1.15 fm coincides roughly with the assumed radius 



of the nucleon of about 1 fm |24j. Our calculation contains 6 radial basis functions, which 
is a sufficient approximation to the solution of Eq.|12|. The convergence of the solution is 
depicted in Fig.0. 

The calculated electric form factor of the proton G P E as a function of the negative photon 
momentum squared is shown in Fig.[|. For A = 0.18 fm we find a very good agreement with 
the experimental data (taken from |2J|). Even for q 2 as high as -3 GeV 2 the calculation fits 
the data well. From the slope of G p E (q 2 ) at q 2 = we obtain for the square root of the mean 
square charge radius 



" 2 ) P 



-6 



dG p E (, 
dq 2 



0.88 fm 



(39) 
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compared to the experimental value y (r 2 ) p = (0.862 ± 0.012) fm pB| . The electric neutron 
form factor is compared to experiment [p7| , p8| in Fig.|5|. We obtain a rather good description 
of the data, with a negative mean square charge radius 

(r 2 ) n = -0.28 fm 2 , (40) 

which is slightly larger than the experimental (r 2 ) n = (— 0.117±0.002) fm 2 0. The magnetic 
form factors of the proton and the neutron are shown in Fig.|| Our results are a factor of 
about 2 too small. However, the calculation matches the experimental data qualitatively. 
The shape of the curves is correct, and we obtain the right signs. Of course, since Gm is 
connected with the spin current, effects of the vector component of a two-quark subsystem 
should be important. Having neglected it, it is even more surprising to find such a good 
correspondance. Looking at Gm{Q), we obtain for the ratio of the magnetic moments of 
neutron and proton 

^ = -0.5 , (41) 

compared to the experimental value —0.685. Of course, the above ratio arises from the 
different charges of the single quark in the nucleons; the scalar diquark does not contribute 
to the magnetic form factor. The prediction of the right ratio jf- is one of the successes of 
the SU(6) quark model. Thus, it would be interesting to study the effects of an additional 
axial- vector diquark component. It should be noted that the neutron form factors are related 
to the differences between the currents 

■n _ -diquark ^ -quark ( A r )\ 

Thus, a description of the neutron form factors is most sensible to the wave function. 

In Fig.0 we present the electric form factor of the proton for three different diquark 
extensions A. The qualitative structure of the curve remains. Fig|] shows the variation of 
the electric form factor for three different values of the diquark mass, with constant A = 0.18 
fm. Obviously, the variation is small. However, we find that the diquark mass should not 
be smaller than about 620 MeV: for ma = 600 MeV the oscillator parameter is a = 1.5 fm, 
and rises sharply with decreasing rrid- 



V. SUMMARY AND OUTLOOK 



We developed a simple, fully relativistic quark-diquark model of the nucleon. In a first 
step, only scalar diquarks are taken into account. Neglecting retardation effects of the 
interaction kernel, a Salpeter equation for a bound state of a constituent quark and a scalar 
particle is obtained. We assume a quark exchange as the only interaction. The electric 
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and magnetic form factors of the nucleon are then calculated in a straightforward way 
using the Mandelstam formalism. With four parameters (constituent quark and diquark 
masses, diquark extension A, diquark-quark coupling g dqq , see Tab.|), we obtain a very good 
description of the electric form factor up to q 2 = — 3 (GeV/c) 2 . The neutron electric form 
factor is reproduced almost quantitatively. The magnetic form factors are in qualitative 
agreement with experiment. 

Obviously, a relativistic treatment is extremely important to describe observables of the 
nucleon in the medium energy range. Furthermore, the assumption of a diquark structure 
of the nucleon seems to be not out of place. 

Of course, it would be interesting to study the effects of eventual axial-vector diquarks. 
Especially for the magnetic form factors they should be important. - The model will be 
extended to baryons with strangeness. Then, various electromagnetic processes may be 
studied. An example is the photoproduction of kaons. The contributing amplitudes can 
be calculated in the framework of the Mandelstam formalism, and a direct coupling of the 



photon to the internal quark- and diquark lines is assumed (see |3C 
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TABLES 



m g 


m d 


9dqq 


A 


350 MeV/c 2 


650 MeV/c 2 


14.15 


1.8 fm 



TABLE I. The parameters of the model 
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FIGURES 




FIG. 1. The Bethe-Salpeter equation for a bound state of a quark and a scalar diquark with a 
quark exchange interaction 
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FIG. 3. The lowest positive eigenvalue of Ti (Eq.|l2|) as a function of the oscillator parameter 
a for different numbers of basis functions (1 to 6 from top to bottom) 
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FIG. 4. The electric form factor of the proton G p E {q 2 ); experimental data are taken from 
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FIG. 5. The electric form factor of the neutron G^{q ); experimental data are taken from [27,28] 
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FIG. 6. The magnetic form factors G^j(q ) of the proton (upper curve) and of the neutron 



(lower curve); experimental data are taken from f25,27 
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FIG. 7. The electric form factor of the proton G p E (q 2 ) for three different values of the diquark 
parameter A 
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FIG. 8. The electric form factor of the proton G p E (q 2 ) for three different values of the diquark 
mass 
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APPENDIX A: TRANSFORMATION PROPERTIES 



Under a Lorentz transformation A, the Bethe-Salpeter amplitude for a bound state of a 
scalar particle and a fermion transforms like a Dirac spinor 



Xp(p) = s a Xap(^p) ■ 



(Al) 



Using S (Ap) = S A S (p)S A we obtain 



T P (p) = SX 1 T AP (A P ) 
Tp( P )=T AP (A P )S A 



(A2) 



Under time reversal (using the standard convention like e.g. in [29|) we find 

Xp(p) = (~l) S - s (~e) 'xp-'iP) 

with S = \ for the nucleon and e = 27 1 7 3 7°. 
Under parity transformation we have 

Xp(p) = ^p1°Xp(P) , 



(A3) 



(A4) 



and under charge conjugation: 



x b P Lh {p) = Scxp l2 {-pY 



bif 2 



(A5) 



with Sc = ij 2 
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